Geometry of quantum dynamics and a time-energy uncertainty relation for mixed states 
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In this paper we establish important relations between Hamiltonian dynamics and Riemannian structures 
on phase spaces for unitarily evolving finite level quantum systems in mixed states. We show that the energy 
dispersion (i.e. llh times the path integral of the energy uncertainty) of a unitary evolution is bounded 
from below by the length of the evolution curve. Also, we show that for each curve of mixed states there 
is a Hamiltonian for which the curve is a solution to the corresponding von Neumann equation, and the 
energy dispersion equals the curve's length. This allows us to express the distance between two mixed states 
in terms of a measurable quantity, and derive a time-energy uncertainty relation for mixed states. In a final 
section we compare our results with an energy dispersion estimate by Uhlmann. 
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I. INTRODUCTION 

Ever since the advent of general relativity, scientists have 
been looking for geometrical principles underlying phys- 
ical laws. Nowadays it is well known that geometry af- 
fects the physics on all length scales, and physical theory 
building consists to a large extent of geometrical considera- 
tions. This paper concerns geometric quantum mechanics, 
a branch of quantum physics that has received much atten- 
tion lately (which is largely due to the crucial role geome- 
try plays in quantum information and quantum computing 
|T]|2)). Here we equip the phase spaces for unitarily evolv- 
ing finite level quantum systems with natural Riemannian 
structures, and establish remarkable but fundamental rela- 
tions between these and Hamiltonian dynamics. 

A quantum system prepared in a pure state is usually 
modeled on a projective Hilbert space, and if the system is 
closed its state will evolve unitarily in this space. Aharonov 
and Anandan 1 3 1 showed that for unitary evolutions there 
is a geometric quantity which, like Berry's celebrated phase 
IU|5], is independent of the particular Hamiltonian used to 
transport a pure state along a given route. More precisely, 
they showed that the energy dispersion (i.e. 1 Ih times the 
path integral of the energy uncertainty) of an evolving state 
equals the Fubini-Study length of the curve traced out by 
the state. Using this, Aharonov and Anandan gave a new 
geometric interpretation of the time-energy uncertainty re- 
lation. 

The state of an experimentally prepared quantum system 
generally exhibits classical uncertainty, and is most appro- 
priately described as a probabilistic mixture of pure states. 
It is common to represent mixed states by density operators, 
and many metrics on spaces of density operators have been 
developed to capture various physical, mathematical, or in- 
formation theoretical aspects of quantum mechanics |6][7]. 
In this paper we utilize a construction by Montgomery 1 8 
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to provide the spaces of isospectral density operators with 
Riemannian metrics, and we show that these metrics admit 
a generalization of the energy dispersion result of Aharonov 
and Anandan to evolutions of finite dimensional quantum 
systems in mixed states. Indeed, we show that the energy 
dispersion of an evolving mixed state is bounded from be- 
low by the length of the curve traced out by the density oper- 
ator of the state, and we show that every curve of isospectral 
density operators is generated by a Hamiltonian for which 
the energy dispersion equals the curve's length. The latter 
result allows us to express the distance between two mixed 
states in terms of a measurable quantity, and we use it to 
derive a time-energy uncertainty principle for mixed states. 

Uhlmann (9||T0l was among the first to develop a mathe- 
matical framework similar to the one presented here. In (TT| 
he used it to derive an estimate for the energy dispersion of 
an evolving mixed state. We compare our energy dispersion 
estimate with that of Uhlmann in this paper's final section. 



II. GEOMETRY OF ORBITS OF ISOSPECTRAL DENSITY 
OPERATORS 

In this paper we will only be interested in finite dimen- 
sional quantum systems that evolve unitarily. They will be 
modeled on a Hilbert space J€ of unspecified dimension 
n, and their states will be represented by density operators. 
Recall that a density operator is a Hermitian, nonnegative 
operator with unit trace. We write Q>{J€) for the space of 
density operators on J€ . 



A. Riemannian structure on orbits of density operators 

A density operator whose evolution is governed by a von 
Neumann equation remains in a single orbit of the left con- 
jugation action of the unitary group of J€ on < 2}{J€). The 
orbits of this action are in one-to-one correspondence with 
the possible spectra for density operators on J€, where by 
the spectrum of a density operator of rank k we mean the 



2 



decreasing sequence 

0= {pi,P2,---,Pk) (1) 

of its, not necessarily distinct, positive eigenvalues. 
Throughout this paper we fix a, and write Qs{o) for the 
corresponding orbit. 

To furnish @(ct) with a geometry, let L(C fc ,<^f) be the 
space of linear maps from C l to jf equipped with the 
Hilbert-Schmidt Hermitian inner product, and P{a) be the 
diagonal k x k matrix that has a as its diagonal. Inspired by 
Montgomery 8 1 , we set 

5fi[q) = {W e UC k ,je) : ¥ tx F = P{a)}, 

and define 

it-.SPia)^ 2(a), ¥ — W 1 . 

Then n is a principal fiber bundle with right acting gauge 
group 

°l£{o) = {U e^{k) : UP [a) = P{a)U], 
whose Lie algebra is 

u(a) = {f eu(Jt):fP(ff) = P(a)0. 

The real part of the Hilbert-Schmidt product restricts to a 
gauge invariant Riemannian metric G on £f[a), 

G(X, Y) = - Tr(X f Y + F + X), 

and we equip @> [a) with the unique metric g that makes n a 
Riemannian submersion. 



B. Mechanical connection 

The vertical and horizontal bundles over £^{a) are the 
subbundles N£fi(a) = Ker tt» and H5f{a) = VSfi(a) L of the 
tangent bundle of Sf[a). Here n* is the differential of n and 
1 denotes orthogonal complement with respect to G. Vec- 
tors in VSf[o) and H5^ [a) are called vertical and horizontal, 
respectively, and a curve in £f(o) is called horizontal if its 
velocity vectors are horizontal. Recall that for every curve p 
in @(cr) and every in the fiber over p(0) there is a unique 
horizontal lift of p to S/'io) that extends from TV This lift 
and p have the same lengths, since n is a Riemannian sub- 
mersion. 

The infinitesimal generators of the gauge group action 
yield canonical isomorphisms between u(cr) and the fibers 
in V^(cr): 

u((i)3(»f(eV ¥ y((T). (2) 

Furthermore, HS^{a) is the kernel bundle of the gauge in- 
variant mechanical connection form si^ = I^ 1 , where Ivp : 
u(ct) — «■ u(cr)* and : TySP[a) — ► u(cr)* are the moment of 
inertia and moment map, respectively, 



The moment of inertia is of constant bi-invariant type since 
it is an adjoint-invariant form on u(cr) which is independent 
of T 1 in 5"(cr). To be exact, 

n*f-i7=i'&((f t t/ + i7 t f)p(a)). (3) 

Using {3} we can derive an explicit formula for the connec- 
tion form. Indeed, if m\, m2,. ..,mi are the multiplicities of 
the different eigenvalues in a, with m\ being the multiplic- 
ity of the greatest eigenvalue, m2 the multiplicity of the sec- 
ond greatest eigenvalue, etc., and if for j = 1,2,...,/, 

Ej = diag(O mi , . . . , Om^j ,l mj , m . +1 , . . . , m; ), 

then 

j 

-\^[LEjX^Ej(-^Ej^XEj) 
= -Tr(X +, P£- l pP t X) 

for every X in T^£^{a) and every £ in u(cr). Hence 

.safy(X) = Y,Ej^XEjP{(j)~ l . 

i 

Observe that the orthogonal projection of Ty5f{o~) onto 
V<p^(cr) is given by the connection form followed by the in- 
finitesimal generator (2). Therefore, the vertical and hori- 
zontal projections of X in T^^tcr) are X L = T'.s/vptX) and 
X 11 = X - Wsfy (X), respectively. 

in. GEOMETRICAL UNCERTAINTY ESTIMATES 

If A is an observable on J€, the gauge invariant vector 
field X^ on S^{a) is defined by 

X A m = -^[exp(^AW] . 
de L V in / J£=o 

Let Xa be the projection of X^ onto @(er), and define the 
uncertainty of A to be the scalar field A A on @ (cr) given by 

AA(p) = Y/Tr(i 2 p)-Tr(ip) 2 . 

We will show that 

AA(p) > h s/g{X A {p),X A (p)), (4) 
AA(p) = h x /g(X A {p),X A {p)) ifX A m e Hy&ia), (5) 

where *P is any element in the fiber over p. 

Assertion (5) follows immediately from the observations 

TrU 2 p) = ^ 2 G(X 1 m,X 1 (^)), (6) 
Tr(ip) = ihTr{^{X A m)P{a)). (7) 
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For if X A C¥) is horizontal, then the right hand side of (6) 
equals h 2 g(X A m,X A m), and the right hand side of (7) 
vanishes. If, on the other hand, X A C¥) is not horizontal, 
we must estimate the difference between GCX^ OP^-XrPP)) 

and Tr(Ap) 2 . The identity 

GCXjmxJ-TO) = -Tr{^{X A m) 2 P{a)), 

together with (6) and (7) yield 

A,4(p) 2 = /l 2 g(X A (p),X A (p)) 

+ ?z 2 Tr(^(^m)P(cr)) 2 (8) 
-h 2 Tr{^{X A m) 2 P{(T)). 

Now (4) follows from the fact that the difference between the 
last two terms in (8) is nonnegative. To see this let U in ^ (cr) 
be such that i Usd^ (XfiQ¥))lfl is a diagonal matrix, say 

iU^{X A m)U f = diag(Ai, A 2 , .... A fc ). 

Then 

Tr(^(X i (¥))P((j)) 2 = -(^P/A;) 2 

= Tr(^(X A m) 2 P(CT)), 
since U commutes with P{a) and x >-» x 2 is convex. 

A. Distance and energy dispersion 

The distance between two density operators with com- 
mon spectrum cr is defined to be the infimum of the lengths 
of all curves in @(cr) that connects them. We will show that 
for any two density operators p and pi in @(cr), 

dist(p ,pi) = iinf C\m P )At, (9) 
n h Jto 

where the infimum is taken over all Hamiltonians H for 
which the boundary value von Neumann equation 

p = X H (p), pUo) = Po, p(fi) = pi, (10) 
is solvable. 

The length of a curve p in @> (cr), with domain to < t < t\ is 

Length[p] = ( \/g{p,p)dt. 

If p = Xtfip), for some Hamiltonian then, by (4), the 
length of p is a lower bound for the energy dispersion: 

Length[p]< \ [ ' AH{p)dt. (11) 

There is a Hamiltonian H that generates a horizontal lift 
of p, because the unitary group of J€ acts transitively on 



L(C fc ,,^f). For such a Hamiltonian we have equality in (TTJ. 
Moreover, we can take p to be length minimizing, in the 
sense that Length[p] = dist(po,pi), because @(cr) is com- 
pact and hence g is complete. Then, 

dist(p ,pi) = \ [ l Aff(p)df, (12) 
by (5). Assertion (9) follows from fTT) and {12}. 

B. Time-energy uncertainty relation 

Consider a quantum system with Hamiltonian H, and 
suppose p is a solution to (To) . The time-average of AH is 

i r f i 

<AH> = — I AH(p)df, At=t!-t . 
At J t(1 

We will show that if po and p\ are distinguishable fl2l[T3l , 
then 

nh 

(AH)At> —. (13) 

For density operators representing pure states, this reduces 
to the time-energy uncertainty relation in (3). 

Let ¥0 in 7r _1 (po) and ¥1 in ^ _1 (pi) be such that 
dist(po,pi) = distcyo.^i)- The operators po and pi have or- 
thogonal supports, being distinguishable, and the same is 
true for T'o and ¥1 since the the support of ¥0 equals the 
support of po, and likewise for and p\. A compact way to 
express this is 

^^1=0, v F|T'o = 0. (14) 

If we consider T'o and *Pj elements in <5?(C k ,Jif), the unit 
sphere in L(C k ,Jtf), they are a distance of nl2 apart. In fact, 
¥(r) = cos{f) x V + sin(r)T'i, with domain < t < nil, is a 
length minimizing curve from T'o to *Pi in 5?(C Con- 
sequently, 

dist(p ,pi) >jt/2. (15) 

The uncertainty relation (13) follows from fTT] and (15) . 
Also note that the estimate (T5) cannot be improved. Di- 
rect computations using fl4) yield T^rO^Cf) = P(cr) and 
, P(r) t T'(fl = 0. Therefore, T'(r) is a horizontal curve in SP (cr), 
and hence (T5) is, in fact, an equality. 

IV. UHLMANN'S BUNDLE AND THE BURES DISTANCE 

Uhlmann 1 1 proved that for unitarily evolving quantum 
systems represented by invertible density operators, the en- 
ergy dispersion is bounded from below by the Bures dis- 
tance between the initial and final states. This result, to- 
gether with ((9) , shows that on orbits of invertible density op - 
erators, the distance function associated with g is bounded 
from below by the Bures distance. Here we present an in- 
dependent argument for this fact, and we give an example 
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of two isospectral density operators between which the two 
metrics measure different distances. 

Let 5^ nv (C", t ^f) be the space of all invertible maps in 
L(C",^f) with unit norm, and @i nv (,^f) be the space of 
all invertible density operators acting on J€. Then n : 
51 nv (C",^f) - @inv(^f) defined by n(f) = is a °lt{n)- 
bundle, which we call Uhlmann's bundle since it first ap- 
peared in 1 10 1 . The geometry of Uhlmann's bundle has been 
thoroughly investigated, and it is an important tool in quan- 
tum information theory, mainly due to its close relationship 
with the Bures metric [141IT51 . 

Uhlmann's bundle is equipped with the mechanical con- 
nection, which means that the horizontal bundle is the or- 
thogonal complement of the vertical bundle with respect 
to the Hilbert-Schmidt metric. Moreover, the metric on 
@inv(^) obtained by declaring n to be a Riemannian sub- 
mersion, is the Bures metric (15| . We denote the associated 
distance function by dists- 

Suppose k — n in (T). Then Sf[ff) is a submanifold 
of 5*5nv(C",^f). Moreover, the vertical bundle of SP{a) 
is subbundle of the restriction of the vertical bundle of 
5^ nv (C" , J€) to Sf{a). However, no nonzero horizontal vec- 
tor in Uhlmann's bundle is tangential to SP(p). To see this, 
let W be any element in SP{o). Then X in T T 5^ nv (C",^) is 
horizontal, i.e. is annihilated by the mechanical connection 
of the Uhlmann bundle, if and only if (TO) 

x f> i X-X i '¥ = 0. (16) 

On the other hand, every X in TySP{a) satisfies 

x V i X+X fx V = (17) 

since = P(p). Clearly, only the zero vector satisfies 
both (16) and (17) ■ 

The distance between po and pi in @>{a) is never smaller 
than Bures distance between them. Indeed, every curve 
between 7T _1 (po) and jr _1 (pi) in SP{a) is a curve between 
n _1 (po) and n _1 (pi) in Sf^iXJ 1 and since the metrics 
on the total spaces of the two bundles are induced from a 
common ambient metric we can conclude that 

dist(p ,pi) > dist B (po,pi)- (18) 

Uhlmann [15 1 and Dittmann HUE] have derived explicit 
formulas for the Bures metric for density operators on finite 
dimensional Hilbert spaces. For density operators on C 2 the 
formula reads 

dist B (p,p + <5p) 2 = ^Tr{8p8p+-^-(8p-p5pf). (19) 



We use (T9) to show that there are density operators po and 
pi acting on C 2 for which the inequality in (T8) is strict. 
Suppose a - [p\,pz), and let e > 0. For e small enough, 

P2 sin 2 (£r) + pi cos 2 (er) (p2 _ Pi) sin(er) cos(e?) 
(P2 — Pi)sin(ef)cos(ef) pi sin 2 (ef) + p2 cos 2 (et) 

is a length minimizing curve in @(cr) between po = p(0) and 
pi - p(l). Thus 

dist(po.pi) = Length[p] = e. 
However, (T9) yields 

dist B (p ,Pi) = I sineli 2 + {P ] - ^ sin 2 e. 

\/2 V 2 PiP2 



V. CONCLUSION 

In this paper we have utilized a construction due to Mont- 
gomery, to equip the spaces of isospectral density operators 
acting on a finite dimensional Hilbert space with Rieman- 
nian metrics, and we have established important relations 
between these and Hamiltonian quantum dynamics. In- 
deed, we have proved that the energy dispersion of a unitar- 
ily evolving density operator is bounded from below by the 
length of the curve traced out by the operator, and that ev- 
ery curve of isospectral density operators can be generated 
by a Hamiltonian such that the energy dispersion equals the 
curve's length. These facts allowed us to express the dis- 
tance between two density operators in terms of a measur- 
able physical quantity, and the paper culminated in a time- 
energy uncertainty estimate for mixed states. In a final sec- 
tion we have compared our energy dispersion results with 
an energy dispersion estimate by Uhlmann. 

We believe that our results have very interesting applica- 
tions in optimal quantum control. Such aspects of the the- 
ory developed here will be investigated by the authors in a 
forthcoming paper. There we will focus on the geometry of, 
and dynamics in, orbits of invertible density operators, and 
we will classify the Hamiltonians that drive states along evo- 
lution curves with minimal energy dispersion. 
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